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The interleukin-1 (IL-1) gene polymorphisms have been impUcated in chronic obstructive pulmonary 
disease (COPD) risk, but results are controversial. We aimed to conduct a meta-analysis to address this 
issue. Odds ratio (OR) and 95% confidence interval (CI) were used to investigate the strength of the 
association. The meta-analysis revealed no association between the ILIB (—511), (—31), (+3954) 
polymorphisms and COPD risk. However, stratification by ethnicity indicated that the T allele carriers of 
the ILIB (—511) polymorphism and the C allele carriers of the ILIB ("—31^ variant were associated with an 
increased risk for developing COPD in East Asians (OR= 1.61, 95% CI: 1.13-2.31, =0.009 and OR= 1.55, 
95% CI: 1.14-2.1 1, =0.006, respectively). The meta-analysis revealed a significant association between the 
ILIRN (VNTR) polymorphism and COPD risk in all study subjects and East Asians under homozygote 
model (22 vs. LL: OR=3.16, 95% CI: 1.23-8.13, P^=0.017 and OR=3.20, 95% CI: 1.13-9.12, P^=0.029, 
respectively). Our meta-analysis suggests that the ILIB (—511), (—31^ and ILIRN (VNTR) polymorphisms 
are associated with COPD risk in East Asians. There is no association between the ILIB ( +3954) 
polymorphism and COPD risk. Further studies should be performed in other ethnic groups besides East 
Asians. 



Chronic obstructive pulmonary disease (COPD) is a chronic inflammatory disease characterized by the 
gradual progression of irreversible airflow obstruction and increased inflammation in the airways and 
lung parenchymal It is a leading cause of chronic morbidity and mortality worldwide. Although cigarette 
smoking has been showed to be a major environmental risk factor for COPD, only 10-20% of smokers develop 
COPD, implying that apart from environmental features, additional risk factors such as genetic variation con- 
tributes to COPD susceptibility^. Genetic linkage studies, candidate gene association studies and genome-wide 
association studies (GWAS) have identified genes that may have roles in the pathogenesis of COPD, including 
microsomal epoxide hydrolase (EPHXl), glutathione S-transferase (GST), interleukin-6 (IL-6), iron-responsive 
element binding protein 2 (IREB2), matrix metallopeptidase 9 (MMP9) and transforming growth factor- P (TGF- 
j^/'^. A better understanding of the genetic architecture of the disease will be important in helping to unravel the 
pathogenetic mechanism and develop novel therapeutic strategies for COPD. 

IL- 1 is a pro-inflammatory cytokine and key contributor to immune responses. IL- 1 occurs in two forms, IL- 1 A 
and IL-IB, both of which bind to the IL-1 receptor^. The IL-1 receptor antagonist (IL-IRA) does not transmit 
any signal and functions as a cell bound inhibitor to both IL-1 A and IL-IB^. The loci for human ILIA, ILIB, 
and ILIRN gene are found as a cluster on chromosome 2ql2 to 2ql4^. The genes are arranged with ILIA 
situated 5-prime, and then ILIB and finally ILIRN 3-prime. Several common polymorphisms within the IL-1 
gene complex have been described, including single nucleotide polymorphisms (SNPs) in ILIB at position 
-511 (rsl6944) and -31 (rsll43627) in the promoter region and at position +3954 (rsll43634) in exon 
5, and a penta-allelic polymorphism representing a variable number of an 86-bp tandem repeat in intron 
2 (rs2234663) of the ILIRN gene^'^. Many genetic association studies have been conducted to investigate 
the relationship of these variants with COPD risk; however, small sample size and varying population 
characteristics led to conflicting results. To validate the potential association between the ILIB, ILIRN poly- 
morphisms and COPD risk, we performed a meta-analysis of data reported in 12 case-control studies from 11 
publications. 
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487 potentially relevant publications 
identified after search for electronic 
databases 



474 publications excluded after 
initial screen of title and abstract 



13 publications for full-text evaluation 



2 publications were excluded owing 
to overlapped sample 



12 studies from 11 publications for final 
analysis 



Figure 1 | Flowchart of study selection process. 

Methods 

Search strategy and study identification. We conducted electronic searches, not 
limited to the English language, of PubMed, Embase, Scopus, Cochrane database, 
China National Knowledge Infrastructure (CNKI) and Wanfang databases to identify 
relevant studies reporting on the association between the ILIB, ILIRN 
polymorphisms and COPD risk in the medical literature from January 1990 up to the 
end of June 2014. Search terms included "chronic obstructive pulmonary disease", 
"COPD", "interleukin-l", "IL-l", "polymorphism", "genetics", and "association". 
References from retrieved papers were also considered. Association studies were 
included in this meta-analysis if they met the following criteria: 1) studies on human 
subjects; 2) studies written in English or Chinese; 3) case-control design; 4) sufficient 
data for examining an odds ratio (OR) with 95% confidence interval (CI). The major 
exclusion criteria were as follows: 1) studies on animal populations; 2) no usable data 
reported; 3) duplicate data; 4) studies presenting deviation from Hardy- Weinberg 
equilibrium (HWE) in the control group. All relevant publications identified through 
the search were scanned on the basis of title, keywords and abstract, and were rejected 
in the initial screening if the paper clearly did not meet the inclusion criteria. Where a 
title/abstract could not be rejected with certainty, the full text of the publications was 
obtained for evaluation. 

Data extraction. Two independent investigators extracted data from the published 
studies and entered them in a customized database. Disagreement was resolved by the 
evaluation of a third reviewer and discussion until a consensus was reached. For each 
study, the following data were collected: first author's name, country, year of 
publication, ethnicity, total number of cases and controls, characteristics of the 
control group, and genotypic frequencies of the IL-l polymorphisms. 

Quality score assessment. The quality of selected studies was evaluated by scoring 
according to Newcastle Ottawa Scale (NOS) (www.ohri. ca/programs/ 
clinical_epidemiolog)^/oxford.asp). This scale consists of three parts relating to 
selection, comparability and ascertainment of exposure. Quality scores ranged from 0 
to 9 and studies were regarded as "high quality" if the score was > 5. Otherwise, 
studies was regarded as "low quality". 

Data analysis. Data were entered and analyzed with the Stata version 1 1.0. We used 
raw data of genotypic distribution, without adjustment for calculation of the study- 
specific estimates of OR and 95% CI. Dominant, recessive and homozygote models 
were evaluated. Z-test was used for assessing the significance of the pooled OR, with 
P<0.05 considered statistically significant. Between- study heterogeneity was 
evaluated with Cochran's Q-test statistics among the studies, with significance level 
set at 0.10. The DerSimonian Laird random- effects model rather than the Mantel- 
Haenszel fixed-effect model was employed to calculate the pooled OR and 95% CF'^°. 
For the ILIRN (VNTR) polymorphism, L signifies any long allele embracing allele 1, 3, 
4, or 5. We created a funnel plot to assess potential publication bias by plotting natural 
logarithm of individual study effect size against the standard error of the natural 
logarithm of individual study effect size. We also assessed publication bias using the 
Begg's test and Egger's test. The HWE law states that if two alleles, T and C, with 
frequencies p and q=l-p, are in equilibrium in a population, then the proportion of 
people with genotypes TT, TC and CC will be p^, 2pq and q^. HWE in the control 
group for each study was tested by using a web-based programme (http://www.oege. 



org/software/hwe-mr-calc.html). Since we only utilized previously published data, we 
did not obtain approval of an ethics committee or written informed consent. 

Results 

Included studies. Figure 1 showed the process of identifying eligible 
studies. 487 publications were identified from the initial keywords 
search. After review, 474 were excluded. Overall, 12 studies from 11 
publications met our inclusion criteria with a total of 1530 COPD 
patients and 1524 controls^^"^\ Table 1 summarized the characte- 
ristics of the studies included in the meta-analysis. In terms of 
ethnicity, eight studies including 1055 cases and 1199 controls 
were performed in East Asians^ two studies containing 
169 cases and 97 controls were undertaken in Arabians ^^'^^ one 
study with 102 cases and 20 controls was conducted in Europeans^^, 
and one study including 204 cases and 208 controls was performed in 
South Asians^^. The publication of Hegab et al. contained an Egyptian 
study and a Japanese study, respectively^^. It was noteworthy that we 
did not extract data for the ILIB (—511), (—31 ) polymorphisms from 
the Liu et al. study since they further evaluated these two polymor- 
phisms using larger sample numbers in another study^^. In addition, 
the data for the ILIB (-31) SNP from the Liu et al. study'^ and the 
data for the ILIRN (VNTR) polymorphism from the Hegab et al. 
study^^, the Lee et al. study^^, and the Shukla et al. study^^ pre- 
sented deviations from HWE in controls, respectively, so we did 
not include them in the final analysis. 

Individual polymorphism meta-analysis. Table 2 showed geno- 
typic distribution of the four polymorphisms. The genotype-wise 
meta-analytic results for each of the four polymorphisms were 
showed in Table 3. 

For the ILIB (—511) polymorphism, nine studies from eight pub- 
Ucations with 1270 cases and 1220 controls were included in the 
meta-analysis^^"^^'^^'^^'^^"^°. There was no significant association 
between the ILIB (—511) polymorphism and COPD risk in all study 
subjects under dominant model (OR=1.22, 95% CI: 0.84-1.75, 
P;, = 0.002, 0.293) (Table 3 and Fig. 2), recessive model 
(OR=0.94, 95% CI: 0.67-1.32, P/,=0.010, P^=0.730) (Table 3), and 
homozygote model (OR=1.10, 95% CI: 0.65-1.85, P/,<0.001, 
P^= 0.734) (Table 3). In subgroup analysis stratified by ethnicity, 
we found that the ILIB (—511) polymorphism was associated with 
COPD risk in East Asians under dominant model (OR=1.61, 95% 
CI: 1.13-2.31, P/, = 0.061, P^=0.009) (Table 3 and Fig. 2), but not 
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under recessive model (OR=1.13, 95% CI: 0.74-1.73, P;, = 0.017, 
P^=0.574) (Table 3) and homozygote model (OR=1.62, 95% CI: 
0.91-2.88, P/,=0.005, P^=0.105) (Table 3). Since there was only 
one study performed in Arabians, Europeans and South Asians, 
respectively, we did not conduct subgroup analysis in these ethnic 
groups. Between-study heterogeneity for the genotype-wise ORs was 
found in dominant model (P= 0.002) (Table 3), recessive model 
(P=0.010) (Table 3), and homozygote model (P<0.001) (Table 3). 

For the ILIB (—31) polymorphism, four studies from three pub- 
lications with 534 cases and 632 controls were included in the meta- 
analysis^^'^^'^^. Pooling data provided no evidence of a relationship 
between this polymorphism and COPD risk in all study subjects in 
dominant model (OR=1.25, 95% CI: 0.79-1.96, P;, = 0.101, 
P^= 0.340) (Table 3 and Fig. 3), recessive model (OR=0.80, 95% 
CI: 0.52-1.21, Ph=0225, P^=0.287) (Table 3) and homozygote 
model (OR=1.00, 95% CI: 0.52-1.95, P;, = 0.065, P^=0.993) 
(Table 3). However, in subgroup analysis based on ethnicity, we 
found a significant association between this SNP and COPD risk 
in East Asians in dominant model (OR=1.55, 95% CI: 1.14-2.11, 
P/, = 0.842, P^= 0.006) (Table 3 and Fig. 3), but not in recessive model 
(OR=0.87, 95% CI: 0.64-1.18, P;,=0.469, P^=0.356) (Table 3) and 
homozygote model (OR=1.27, 95% CI: 0.87-1.87, P/, = 0.599, 
P^^ 0.221) (Table 3). Between-study heterogeneity was found for 
the genotype-wise OR in homozygote model (P= 0.065) (Table 3). 

For the ILIB (+3954) polymorphism, four studies from three 
publications with 406 cases and 359 controls were included in the 



meta-analysis^^'^^'^^. Available data did not suggest an association 
between this polymorphism and COPD risk in all study subjects in 
dominant model (OR=1.16, 95% CI: 0.78-1.73, P/,=0.830, 
P^= 0.464) (Table 3), recessive model (OR=0.97, 95% CI: 0.32- 
2.88, P;,= 0.486, P^= 0.951) (Table 3) and homozygote model 
(OR=1.02, 95% CI: 0.33-3.12, P;, = 0.480, P^=0.974) (Table 3). In 
subgroup analysis according to ethnicity, we also did not find an 
association between this SNP and COPD risk in East Asians 
(Table 3). Because of limited availability of published results, we were 
unable to perform subgroup analysis in Arabians, Europeans and 
South Asians, respectively. Between-study heterogeneity for the 
genotype-wise ORs was not found (Table 3). 

For the ILIRN (VNTR) polymorphism, five studies with 322 cases 
and 362 controls were included^^'^^'^^'^^'^\ The pooled effect estimates 
among all studies suggested a significant association between the 
ILIRN (VNTR) polymorphism and COPD risk in recessive model 
(OR=2.59, 95% CI: 1.02-6.58, P;,=0.891, P^=0.046) (Table 3) and 
homozygote model (OR=3.16, 95% CI: 1.23-8.13, P;, = 0.800, 
P^^ 0.0 17) (Table 3 and Fig. 4), but not in dominant model 
(OR=1.64, 95% CI: 0.99-2.73, P;, = 0.241, P^=0.056) (Table 3). In 
subgroup analysis stratified by ethnicity, we found that the ILIRN 
(VNTR) polymorphism was associated with COPD risk in East 
Asians in homozygote model (OR=3.20, 95% CI: 1.13-9.12, 
P;, = 0.605, P^=0.029) (Table 3 and Fig. 4), but not in dominant 
model (OR=1.48, 95% CI: 0.77-2.83, P;, = 0.146, P^=0.241) 
(Table 3) and recessive model (OR=2.60, 95% CI: 0.93-7.31, 
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Begg's test and Egger's test suggested no evidence of publication 
bias (P= 1.000 and P=0.865, respectively). 

Discussion 

The ILIB and ILIRN polymorphisms have been implicated in the 
pathogenesis of COPD and have been investigated in numerous case- 
control association studies. However, the results are controversial, 
possibly because single studies may have been underpowered. Given 
the amount of accumulated data and the still equivocal role of the 
IL-1 gene polymorphisms in the etiology of COPD, we performed 
the present meta-analysis of 12 case-control studies for a total of 
1530 COPD patients and 1524 controls to investigate the association 
between the four most studied ILIB and ILIRN polymorphisms 
and COPD. The main findings of the present meta-analysis are: 
(1) the T allele carriers of the ILIB (—511) polymorphism are 
associated with an increased risk for developing COPD in East 
Asians; (2) the C allele carriers of the ILIB (—31) polymorphism 
are associated with an increased risk for COPD in East Asians; (3) 
compared with the LL homozygotes, the 22 homozygotes of the 
ILIRN (VNTR) polymorphism are associated with an increased risk 
for developing COPD in East Asians; and (4) there is no significant 
association between the ILIB (+3954) polymorphism and COPD 
risk. 

IL-IB is a pro -inflammatory cytokine which is synthesized by a 
variety of cell types, including blood monocytes and tissue macro- 
phages. IL-IB is thought to be an important mediator in cigarette 
smoke-induced inflammation and COPD. In BALB/c mice, IL-IB 
production was significantly up-regulated in lung homogenates after 
smoke exposure^^. In human subjects, elevated IL-IB levels were 
observed among smokers^^ and in primary explant cultures of bron- 
chial epithelial cells derived from COPD patients^^. IL-IB expression 
in COPD neutrophils was correlated with disease severity as mea- 
sured by forced expiratory volume in 1 second (FEVl)^^. In addition, 
sputum IL-IB was showed to be a potential biomarker for bacteria- 
associated exacerbations of COPD^^. IL-IB plays a central role in the 
regulation of immune responses and inflammatory processes, 
including promotion of the movement of inflammatory cells from 
the blood to inflamed tissues, regulation of the extracellular matrix, 
induction of the expression of a variety of inflammatory mediators 
such as IL-6, IL-8 and tumor necrosis factor- a(TNF- a), and pro- 
motion of the differentiation of inflammatory cells^. Since IL-IB 
regulates a variety of inflammatory processes, any alteration in its 
blood or tissue concentration may significantly affect these processes. 
The ILIB (—511)T and (—31)C alleles are associated with higher 



Table 3 | Meta-analysis of the relationship of the IL-1 polymorphisms with COPD risk 
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CI, confidence interval; COPD, chronic obstructive pulmonary disease; IL-1 , interleukin-1 ; NA, not available; OR, odds ratio; P|^et, P-value for heterogeneity; P^, P-value for overall effect. 
"For ILIB (-51 1) polymorphism: dominant (TT + TC vs CC), recessive (TT vs TC + CC) and homozygote (TT vs CC). 
^For ILTB (-31) polymorphism: dominant (CC + CT vs TT), recessive (CC vs CT + TT) and homozygote (CC vs TT). 
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'^For ILIRN (VNTR) polymorphism: dominant (22 + 2Lvs LL), recessive (22 vs 2L + LL) and homozygote (22 vs LL). 
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Ph=0.732, Pz=0.069) (Table 3). We did not perform subgroup ana- 
lysis in Arabians, Europeans and South Asians because of limited 
availability of published results. Between-study heterogeneity for the 
genotype- wise OR was not indicated (Table 3). 

Publication bias. Since publication bias was hard to detect when the 
number of studies were small, we selected the ILIB (—511) 
polymorphism to assess publication bias (nine studies included). 
The shape of the funnel plot seemed symmetrical (Fig. 5). Both 
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Figure 2 | Meta-analysis for the association between the ILIB (-511) polymorphism and COPD risk in dominant model. Each study is shown by the 
point estimate of the odds ratio, and a horizontal line denotes the 95% confidence interval. The pooled odds ratio is represented by a diamond. The 
area of the grey squares reflects the weight of the study in the meta-analysis. 



levels of IL-IB and with severe inflammation, in comparison to 
(—511)C and (—31)T alleles, which are associated with lower levels 
of IL-IB^^'^^. Genetic variation in the level of gene expression within 
the ILIB gene may alter susceptibility to COPD. In the present meta- 
analysis, we found that the T allele carriers of the ILIB (—511) SNP 
and the C allele carriers of the ILIB (—31) polymorphism had an 



increased risk for COPD in East Asians, suggesting an effect of the 
ILIB gene on COPD in this ethnic group. Due to limited availability 
of published data, subgroup analysis was only conducted in East 
Asians. To test whether this association is population dependent, 
future studies should be performed in other ethnic groups besides 
East Asians. 
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Figure 3 | Meta-analysis for the association between the ILIB ( -31 ) polymorphism and COPD risk in dominant model. Each study is shown by the point 
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area of the grey squares reflects the weight of the study in the meta-analysis. 
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Figure 4 | Meta-analysis for the association between the ILIRN (VNTR) polymorphism and COPD risk in homozygote model. Each study is 
shown by the point estimate of the odds ratio, and a horizontal line denotes the 95% confidence interval. The pooled odds ratio is represented by a 
diamond. The area of the grey squares reflects the weight of the study in the meta-analysis. 



IL-RA binds to the IL-1 receptor and acts as a competitive inhib- 
itor of IL-IB. Allele 2 (two repeats) of the ILIRN (VNTR) poly- 
morphism is associated with increased IL-lra levels^^. In addition, 
allele 2 strongly enhances in vitro production of IL-1 B and plays a 
pivotal role in determining the balance between secreted IL-IB and 
IL-1RA'°. The ILl-RA/IL-lB ratio is critical in determining the 
severity of inflammatory responses. The ILIRN allele 2 has been 
showed to be associated with a low ILl-RA/IL-lB ratio, thereby 
leading to a more prolonged and more severe pro -inflammatory 
immune response^\The frequency of this allele was found to be 
increased in diseases of autoimmune or inflammatory nature, 
including ankylosing spondylitis^^, diabetic nephropathy^^, and 
ulcerative colitis^^. Our meta-analysis showed that the ILIRN allele 



2 homozygotes had an increased risk for COPD, suggesting a signifi- 
cant pathogenetic effect of the ILIRN aflele 2 in COPD. 

There are several limitations to this meta-analysis. First, between- 
study heterogeneity was found in some pooled analyses. We were 
unable to further identify the exact sources of heterogeneity by meta- 
regression analysis in that there were limited relevant data provided 
by the included studies. Some potentially relevant factors, such as 
ethnicity, genotyping method, gender, phenotype of the disease may 
account for heterogeneity. We further performed subgroup analysis 
according to ethnicity, which greatly reduced heterogeneity. Second, 
we only included publications written in English or Chinese in the 
present meta-analysis, which may have missed potentially eligible 
studies published in other languages. Third, few included studies 



Begg's funnel plot \'^/ith pseudo 95% confidence limits 



^ 0 



-1 



s.e. of: logor 

Figure 5 | Begg's funnel plot for the ILIB (-511) polymorphism and COPD risk. 
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reported the percentage of smokers in COPD patients and controls. 
Because of limited availability of published data, we were unable to 
evaluate the effect of smoking on the association between the ILl 
gene polymorphisms and COPD. Fourth, there has been consid- 
erable strength of linkage disequilibrium (LD) in the IL-1 gene clus- 
ter. Feakes et al. found by family studies significant LD between the 
ILIB (-511) and ILIRN (VNTR) polymorphisms'^ Joos et al. found 
that the ILIB (-511)1: allele in combination of the ILIRN allele 2 was 
protective against a rapid decline in lung function in smokers'^. 
These results suggest that haplotype analysis, as well as genotype 
analysis, may be necessary in evaluating genetic effects of the IL-1 
gene cluster on COPD. We did not address this issue in our meta- 
analysis since few included studies performed haplotype analysis and 
there were discrepancies in constructing haplotypes among 
them^^'^^'^^. Association studies may provide more insights for the 
role of the IL-1 gene in COPD if researchers can investigate haplo- 
type association. 

In summary, our meta-analysis showed that the ILIB (—511), 
(—31) and ILIRN (VNTR) polymorphisms are associated with 
COPD risk in East Asians. There is no significant association 
between the ILIB (+3954) polymorphism and COPD risk. Further 
studies with large sample size are needed to establish a more defin- 
itive conclusion. In addition, more studies should be performed in 
other ethnic groups besides East Asians. 
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